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Shell-model calculations for several mass 16-19 nuclei are performed using the N'^LO two-nucleon 
potential V2N with and without the addition of an in-medium three-nucleon potential V-^'^ , which 
is a density-dependent effective two-nucleon potential recently derived from the leading-order chiral 
three-nucleon force Vsjv by Holt, Kaiser, and Weise. We first calculate the Viom-k low-momentum 
interactions from V2N and V^"^ . The shell- model effective interactions for both the sd one-shell 
and sdpf two-shell model spaces are then obtained from these low-momentum interactions using 
respectively the Lee-Suzuki and the recently developed Okamoto and Suzuki iteration methods. The 
effects of V^''' to the low-lying states of ^*0, ^^F, ^^O and ^^F are generally small and attractive, 
mainly lowering the ground-state energies of these nuclei and making them in better agreements 
with experiments than those calculated with V2N alone. The excitation spectra of these nuclei are 
not significantly affected by V^"^ . The low- lying spectra of these nuclei calculated with the sd and 
sdpf model spaces are closely similar to each other. Our shell-model calculations for ^^O indicate 
that the V^"^ interaction is important and desirable for the binding energy of this nucleus. 
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I. INTRODUCTION 

The nuclear shell model has been very successful in 
microscopic descriptions of nuclear structure, and in this 
approach the shell- model effective interaction Veff has 
played an important role and its determination has been 
extensively studied As discussed in these refer- 

ences, Veff may be determined using either an empiri- 
cal approach where it is required to reproduce selected 
experimantal data or a microscopic one where Veff is 
derived from realistic nucleon-nucleon interactions using 
many-body methods. The interactions used in such mi- 
croscopic calculations have been mostly the two-nucleon 
(NN) interaction V2N 0, 0] • Should 14/ / have also con- 
tributions from the three-nucleon (NNN) force V3N in 
addition to the two-nucleon one? In fact the need of V3N 
in nuclear many-body problems has long been recognized. 
The use of V2N alone has been inadequate in reproducing 
the empirical nuclear matter saturation properties (see 
e.g. and references quoted therein). The inclusion 

of VsN has been of essential importance in describing the 
binding energies and low-lying spectra of light nuclie 0- 
Q and in explaining the long half-life of the ^''C — J*^"* N 
/3-decay [13, Hi- Otsuka et al. [l^ have shown that 



the inclusion of V^n plays a crucial role in describing the 
oxygen isotopes near the drip line. 

In the present work, we shall calculate the shell-model 
effective interactions for the sd and sdpf shells using the 
chiral N^LO two-nucleon potential V2N [13| with and 



without the inclusion of the in- medium NNN force Vg'J^'^^, 
which is a density-dependent two-nucleon potential re- 
cently derived from the leading-order chiral NNN force 
V3JV by Holt, Kaiser and Weise 13, 111- We shall apply 



these effective interactions to shell-model calculations for 
nuclei ^^O, ^^F, ^^O and ^^F, to study the effects of 
We shall first calculate the low- 
iM3 from V2N and 



^med these nuclei, 
momentum Viow-k interactions 



Our shell model effective interactions Veff will 
then be calculated from these Viow-k interactions using 



T/med 



a folded-diagram formalism 

@, i, im, ll8| . In this formal- 
ism, Veff is given as a folded-diagram expansion. For 
the degenerate sd case, this expansion can be summed 
up using the commonly employed Lee-Suzuki (LS) it- 
eration method jlj;, 2^- For the non-degenerate sdpf 
case we shall sum up the expansion using the extended 
Krenciglowa-Kuo iteration method recently developed by 
Okamoto, Suzuki, Kumagai and Fujii (EKKO) ^21, 23]. 
The EKKO method is efficient for deriving the effective 
interactions of non-degenerate model spaces such as the 
sdpf two-shell space, and it can also be conveniently ap- 
plied to degenerate model spaces such as the sd one-shell 
one. 
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As mentioned earlier, the in-medium NNN potential 
^med (fependent on the nuclear medium density p. We 
shall study this p dependence for both low and moder- 
ately high densities. Clearly, at zero density V^''' van- 
ishes. Nuclear matter calculations using V2N alone have 
not been able to satisfactorily reproduce the empirical 
nuclear matter saturation properties [3, It may be 
useful to study if the p dependence of V^"^ near the nu- 
clear matter saturation density po (0.16 fm~^) may play 
an important role for nuclear matter saturation. We shall 
do so by carrying out nuclear matter calculations with 
and without the inclusion of V^'^. Similarly this de- 
pendence may also be important for closed-shell nuclei, 
such as ^^O, whose nucleons are embedded in a nuclear 
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medium of densities near po- We shall calculate closed- nucleon i and 
shell nucleus ^^O with and without the inclusion of V^^^^ 
as a further study of its effect near po- The valence nu- 
cleons, such as those of ^®0, are in a medium of densities 
much less than po- Thus our shell- model calculations for 
valence nuclei mentioned earlier are a study of V^'^ at 
low densities. 

The organization of the present paper is as follows. We 
shall describe some details about the derivation of V^'^ 
from V^M in section II. It may be noted that V^'^ is an ef- 
fective density-dependent two-nucleon interaction, which 
is more convenient than its underlying three-nucleon po- 
tential for nuclear many-body calculations. The meth- 
ods we shall employ for the derivation of the shell-model 
effective interactions from V2N and V^''' will be 

outlined there. In section III we shall present first our 
results of nuclear matter calculations using V2N with and 
without the inclusion of VJj'JJf''. The ring-diagram formal- 
ism 0| employed for our nuclear matter calculations 
will also be outlined. Next we shall report the results 
of a similar ring-diagram calculation for ^^O. Results of 
our shell-model calculations of ^^O, ^^F, ^^O and ^^F 
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The low-energy constants ci = — 0.76GeV ^,03 = 
-4.78GeV~\ and C4 = 3.96GeV^ appear already in 
the N^LO two-nucleon potential and are therefore con- 
strained by low-energy NN phase shifts [27]. The low- 
energy constants cd and ce are typically fit to reproduce 
the properties of light nuclei [7|-|9| . 

A general three-body force may be written in second 
quantization as 



3^ = ^ E V{[l2i],[A^Q])b\blblKhzbi, 



(5) 



123456 



where the antisymmetrized matrix element is 



using V2N with and without the inclusion of will 
be presented and discussed in this section. A summary 
and conclusion is presented in section IV. 



II. FORMALISM 

We first describe how we include the effects of the 
leading-order chiral three-nucleon interaction, V^m, in 
our calculations. We consider the nuclear interaction V 
as given hj V ^ {V2N + ^3"^^°!), where V2N is the N^LO 
Idaho two-nucleon potential [13| and V^'^ is a density- 
dependent two-body interaction obtained from the chiral 
three-nucleon force by closing one pair of external lines 
and summing over the filled Fermi sea of nucleons. The 
leading contribution to V-^m occurs at N^LO in the chiral 
power counting and is composed of a long-range two-pion 
exchange component V^^, a medium-range one-pion ex- 
change term V^^, and a pure contact interaction V^l^: 



y([123], [456]) EE (123|V3Arl456+645+564-654-546--465). 

(6) 

Here (123]V3Ar]456) is a simple product matrix element, 
and 6^ and b are creation and destruction operators de- 
fined with respect to the particle- hole vacuum \C) with 
bk\C) = for all k. From eqs. (Il]|4|) we can write 



V'iN — V^pj + V^^ + V^"^ , where V^'j^ is the component 
of VsAT that is symmetric with respect to the interchange 
i k. Now we contract one pair of the and b oper- 
ators of the above V^n (both operators must be holes), 
and this leads to an effective two-body force 



.(2) 



.(3) 



T/mcd 



i^i?([12],[45])6l465&4, 



(7) 
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with 

^([12], [45] 
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(1^21^5^^^ |4i5) 



12^1^3(^)145*)- (1*21^3(^)1*45) 

3;i4*5)- (12*11/3(^)14*5) 
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where 4 5 denotes the nine exchange terms. The above 
result is unchanged when V^fj is replaced by either v!^]^ 

(2) or V^^ . In our calculations we consider a background 
medium of symmetric nuclear matter at constant den- 
sity characterized by a Fermi momentum kp- In this 



way analytic expressions can be obtained for 



3N 



(3) 



where gA = 1.29,/^ = 92.4MeV, A^^ = 700MeV,m^ = 
138.04MeV/c^ is the average pion mass, qi = p/ — pi is 
the difference between the final and initial momentum of 



mcd 

3N ' 

shown in refs. [10|, The above is a density depen- 
dent effective 'two-nucleon' interaction which, unlike its 
underlying three-nucleon force, can be readily used in 
many-body problems. As detailed in flT], the partial- 
wave potentials of the above V^''' have been derived 
from the lowest-order three-nucleon force. We shall use 
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them in our calculations, namely we shall consider the 
nucleon interaction as given by {V2N + YIk^- 
We use the folded-diagram theory [i~ 



17, 181 to cal- 



culate the effective interaction. Briefly speaking, in this 
theory Veff is given by a folded-diagram series 



Veff = Q-Q J Q+Q J Q J Q-Q J Q J Q J Q --, 

(9) 

where Q represents a so-called Q-box consisted of irre- 
ducible diagrams, as illustrated by the 1st- and 2nd-order 
diagrams of Fig. 1. (The Q'-box is the same as the Q-box 
except that Q' does not have diagrams Ist-order in the 
interaction.) Suppose we include only V2n- Then each 
vertex in the diagrams shown represents a Viow-k low- 
momentum interaction [li-16| derived from V2n- When 
we include also V^^'^ in the calculation of Veff, the Q- 
box diagrams for effective interactions will have a specific 
type of V^'"^ vertices as illustrated in Fig. 2. Recall that 
diagram d4 of Fig. 1 represents the interaction between 
two valence nucleons via V2n- This diagram becomes 
diagram (a) of Fig. 2 if its V2N is replaced by V^'^. 
Here the two valence nucleons have to have the partic- 
ipation of a sea-nucleon (below Fermi sea) in order to 
activate V3N, as indicated by the hole-line loop in the 
diagram. ( Note each V^'^'^ vertex of Fig. 2 represents 
a Viow~k interaction derived from V^jlf^.) Similarly dia- 
gram dl of Fig. 1 represents the interaction of a valence 
nucleon with a sea-nucleon via the two-nucleon interac- 
tion V2N- If this interaction is replaced by V^"^, this 
diagram becomes diagram (b) of Fig. 2. V^n must in- 
volve three nucleons, and hence here the valence nucleon 
interacts with 'two' sea-nucleons as indicated by the two 
hole-line loops attached to V^n- Diagram (c) is the V^n 
core polarization diagram corresponding to diagram d7 
of Fig. 1. Again here one hole- line loop is needed for 
each vertex to activate V37V As seen from Eqs. (31) and 
(32), Vg'JJf' requires the involvement of at least one sea- 
nucleon, and it is this requirement which is reflected by 
the hole-line loops attached to the V^n vertices in Fig. 2. 
To summarize, the Q-box diagrams of Fig. 1 are used for 
calculating Veff when we use only V2n, while additional 
diagrams as illustrated by Fig. 2 are also included when 
the interaction is consisted of both V2N and V^n- 

After obtaining the Q-box, we calcualte the effective 
interaction of Eg. (9) using iteration methods. The 
LS method 

Urn 

is convenient for degenerate model 
spaces, and it is used to calculate the sd one-shell effec- 
tive interactions. For the non-degenerate sdpf two-shell 
effective interactions, one can use either the Krenciglowa- 
Kuo iteration method (KK) [2^ 



2J| or the recent ex- 



tended KK method (EKKO) of Okamoto et al. [21 1. 
Both the EKKO and KK methods are convenient for 
calculating the effective interaction for non-degenerate 
model spaces, while the EKKO method being more effi- 
cient (faster converging rate) . 25] . The main difference 
between the EKKO and KK methods is the following. 





FIG. 1: Low-order diagrams constituting the Q-box. 
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FIG. 2: Q-box diagrams with three-nucleon forces. Each 
Viitf'' interaction is represented by a 'narrow band with hole- 
line-loop' vertex. 



The KK method employs the irreducible vertex function 
Q-box which may be written as 



Q{uj) = [PVP + PVQ 



1 



a; - QHQ 



QVP] 



(10) 



where V represents the NN interaction, P denotes the 
model-space projection operator, and Q = {1 — P). The 
energy variable to is determined self-consistently 1^, 2^ . 
Note that Q contains only valence-linked diagrams as in- 
dicated by the subscript L. In the EKKO method, the 
above Q-box is replaced by the Z-box defined by 21 , 2^ 



Ziuj) 



1 



i-Qi(w) 



[QiLj)^QliLu)PiLU~Ho)P], (11) 



where Qi is the first derivative dQ/du. When the P- and 
Q-space are not sufficiently separated, the above Q-box 
may have singularities. An important advantage of the 
EKKO method is that the Z-hox is well behaved even 
when the corresponding Q-box is singular, overcoming 
the above singularity difficulty and making the EKKO 
method more efficient [lHli^. We shall calculate the sdpf 
two-shell effective interactions using both the EKKO and 
KK iteration methods, to cross check their convergences. 
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III. RESULTS AND DISCUSSION 

As described in section II, we shall calculate the shell- 
model effective interactions with the inclusion of the 
medium-dependent three-nucleon force V^'^'^ which is ob- 
tained from a chiral NNN force by integrating one partic- 
ipating sea-nucleon over the Fermi sea. The V^'"^ inter- 
action is a density dependent effective interaction, and 
we shall first carry out nuclear matter calculations with 
the inclusion of V^^lf^, to study the properties of this 
interaction at various densities near the nuclear matter 
saturation density po- Before doing so, we need to choose 
or decide the low-energy constant of V^n (see Eqs. (1- 
3)) to be used in our calculations. The low-energy con- 
stants ci, C3 and C4 of V^^"^ are well determined; they 
are constrained by low-energy NN scattering data [27| . 
Their values so determined (section II) will be used in 
the present work. But the low-energy constants cd and 
cb, of V^*'^'^'' and V^'^'' respectively, are less well known; 
their values determined from properties of light nuclei 
exhibit considerable variations [7H10|. The V^^'^ inter- 
action depends explicitly on the Fermi momentum kp 
which is well defined for nuclear matter. But kp is not 
well defined for finite nuclei, causing uncertainty in deter- 
mining the low-energy constants of VJ^^"^ from properties 
of finite nuclei. But for nuclear matter, kp is well defined; 
we have thus chosen to study these constants by way of 
nuclear matter calculations with V^"^ included. 




FIG. 3: Nuclear matter ring diagrams with vertices from V2N 
(wavy line) and V^'^ (narrow band with hole-line- loop). 

We calculate the equation of state (EOS) for symmet- 
ric nuclear matter using a ring-digram formalism 0, Q . 
A brief description of this formalism is presented below, 
to outline how we include V^^'^ in our calculations. Us- 



ing familiar renormlization procedures l3- 16j|. we first 
calculate the low-momentum interactions V^o^_fe(A) and 
^ilw-kW respectively from V2N [11 and V^Jf^. A com- 
mon decimation scale of A = 2.1fm~^ is employed. This 
value is chosen because at this scale the low-momentum 
interactions derived from different NN potentials [13, 29- 
31 1 are remarkably close to each other [l6|, leading to 
a nearly unique low-momentum interacton. Using the 
above ring-diagram framework, the ground-state energy 
shift AE is given by the all-order sum of the pphh ring 
diagrams as illustrated in Fig. 3. {AEq is defined as 



^0 



and 4'^"'= 



where Eq is the true ground state energy 
that for the non-interacting system.) As shown 



in Fig. 3, diagram (1) is a Ist-order ring diagram with its 
"'^ calculated from the N^LO chiral NN po- 
Similarly diagram (2) is a Ist-order ring dia- 



vertex K- 
tential [1 

gram with V^'^^'^ obtained from the lowest order chiral 



V3N as described in section II. Diagram (3) has vertices 
from both V^'^''' and V^'^'''. Similar to Fig. 2, to have 
the pair of nucleons in the ring diagrams interact with 
VsN there must be the participation of a third nucleon. 
Thus the V^n vertices in Fig. 3 all have a one-hole-line 
loop attached to them. 

With these ring diagram summed to all orders, the 
ground-state energy shift for nuclear matter is given as 

AEo = fdXj2 E Y,n{ij,X)Y:,{kl,X) 

•'^ m ijkKA 

x{imi:^_,iA)+v^i^_,{A)]\ki). (12) 



The transition amplitudes above are Yj^{kl,X) — 
(*,„(A, A - 2)|/3(fcOI*o(A, A)), where *o(A,A) denotes 
the true ground state of nuclear matter which has A par- 
ticles, ^m(A, A — 2) the mth true eigenstate of the {A-2) 
system, and f3{kl) — bibk if {k,l) > kp and = if 
{k,l) < kp. Note that A is a strength parameter, in- 
tegrated from to 1. The amplitudes Y are calculated 
from a RPA equation dSlii] based on the (F2JV+K5T'') 
interactions. Note that our calculation is the same as 
the usual Hartree-Fock (HF) one if we include only the 
Ist-ordcr ring diagrams (1) and (2) of Fig. 4, and in this 
case the above Y amplitude becomes Y{ij) = niJij where 
ni=l for i < kp and =0 otherwise. Our ring-diagram cal- 
culations include particle-hole fluctuations of the Fermi 
sea, while not so for the HF calculations. 
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FIG. 4: Ring-diagram equation of state for symmetric nuclear 
matter calculated with (w/ 3NF) and without (no 3NF) V^^"^. 

We have carried out ring-diagram calculations for sym- 
metric nuclear matter using a wide range of values for cd 
and cp- Our results with cd = —2.7 and cp = 0.7 are 
displayed in Fig. 4; they give Eq/A ~ — 16MeV and 
Pq ~ 0.16fm~^, both in satisfactory agreement with the 
empirical values. Two curves are shown in the figure. 
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one with V2N alone and the other with the addition of 
V^'^. Comparing with the V2N curve, it is of inter- 



est that the effect of Vg'^f 



is sHghtly attractive for low 
densities (<^ 2po/3) while becomes strongly repulsive 
at high densities. Recall that we have used a common 
decimation scale of A = 2.1fm-^ ^ for both V^'^^'' 
and V'g^'""'^. It is of interest that in the ring-diagram 
calculations with Brown-Rho scaling 0, Q a large A of 
~ 3.0fm^^ is needed for obtaining satisfactory nuclear 
matter saturation properties, while in the present cal- 
culation with Kj'JJf'^ satisfactory results can be obtained 
with the use of a smaller A (2.1 fm~^). 

We have calculated the HF s.p. spectrum ek in nuclear 
matter with the inclusion of both V2N and three-nucleon 
force Vg'JJf'' . The inclusion of the latter has been found to 
have significant effect to the spectrum, ek can be well fit- 
ted by the quadratic expression {k'^h^ /{2m*) + A) where 
TO* is the effective mass and A is a well-depth parameter 
representing the s.p. energy at zero momentum. In Fig. 
5 we present our results for A and to/to* for various den- 
sities. Two curves are shown; the lower one without and 
the top one with the inclusion of V^'jlf'^. The densities 
for the 7 data points of each curve are, from left to right, 
(0.25, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 po) respectively. As shown 
by the lower curve of the figure, A and m/m* both vary 
monotonically with the density when only V2N is em- 
ployed. With increasing density, A becomes incresingly 
more negative and m/m* increasingly larger, exhibiting 
no saturation (up to 3.0 po)- The trend is quite differ- 
ent for the upper curve where the three-nucleon force is 
included. Here m/m* still monotonically increases with 
density, but A arises after saturation, approaching zero 
at some higher density. It is of interest that Vg'JJ/^'^ has a 
large effect in raising the chemical potential of nucleons 
in high density nuclear matter; this raise of the chemical 
potential may play an important role in enabling such 
nucleons decaying into other baryons such as hyperons. 
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FIG. 5: Comparison of m* and A for the in-medium s.p. 
spectrum calculated with (w/ 3NF) and without (no 3NF) 



See text for other explanations. 



As indicated by Fig. 4, our ring-diagram calculations 
with the inclusion of the NNN force VJj'JJf' have described 
well the nuclear matter saturation properties. It would be 
useful and of interest to check if this approach is also sat- 
isfactory for describing closed-shell nuclei whose nucleons 
are mainly embedded in a similar nuclear medium of den- 
sities near pQ, the nuclear matter satturation density. To 
study this possiblity, we have extended the above ring- 
diagram nuclear matter calculation to the closed shell 
nucleus ^^O, using the same low-energy parameters Co 
and Ce (-2.70 and 0.70) determined earlier. Note that 
these parametrs will also be used later in our sd and 
sdpf shell-model calculations. Our ring-diagram calcu- 
lation for ^^O is quite similar to the nuclear matter one. 
The ground-state energy shift for ^^O is also given by the 
all-order sum of the ring diagrams illustrated in Fig. 3, 
except for the following differences. For nuclear matter, 
the s.p. states are plane-wave states and each particle 
line in the figure denotes a plane- wave particle with mo- 
mentum k > kp and each hole line with k < kp- For ^^O 
we use harmonic oscillator s.p. wave functions {huj—l-i 
MeV), and each particle line in the figure denotes a nu- 
cleon in the particle shells {IsOdlpOf), and each hole line 
a nucleon in the OsOp shells. For nuclear matter, we use 
momemtum model space where all particles have momen- 
tum fc < A. For ^^O, we use an oscillator model space 
composed of the s.p. orbits of the OsOplsOdlpOf shells. 
Indeed the ring-diagram calculations for nuclear matter 
and ^^O are quite similar. The ground-state energy shift 
for ^^O is also given by Eq.(12), except that the orbits 
{ijkl) now all refer to shell-model orbits and the summa- 
tion restriction of {ijkl < A) is replaced by that of ijkl 
belonging to the above oscillator model space. 

With the same V2N and Vg'JJf'^ interactions used for 
our nuclear matter calculations, we have carried out ring- 
diagram calculations for the ground-state energy of ^^O. 
Since V^'^ is dependent on the nuclear density p, we 
need to have the value of p to carry on the calculation. 
For nuclear matter, this p is well defined, but not so for 
finite nuclei. A preliminary estimate for the density for 
^^O can be made by considering it as a uniform sphere 
[3^ ; in this way the average density for ^^O is estimated 
as ~ 0.85/OO- A more realistic estimate for this density 
should be lower, as physically the nuclear surface is dif- 
fused (not sharp). In Table I we present results for the 
ground-state energy per nucleon of ^^O for sevearl den- 
sities, p/po = 0, 0.5, 0.70 and 0.75. As seen, the binding 
energy calculated with p — (i.e. no V^'^) is clearly not 
adequate. The inclusion of V^^'/f'^ lowers the ground-state 
energy per nucleon significantly; for example the inclu- 
sion of V^^'^ at p/ po=0.75 has increased the ring-diagram 
ground-state energy from -6.97 (only V2n) to -8.01 MeV 
(with V^'^) per nucleon, in good agreement with exper- 
iment. Since we have not included the Coulomb inter- 
action in our calculations, an empirical Coulomb energy 
of 1.14 MeV per nucleon [s^ has been included in the 
above calculation for the ground-state energy of ^^O. It 



6 



TABLE I: Ring-diagram calculations of the ground-state en- 
ergy per nucleon Eq for ^^O. The interaction (V2JV + Vsn"') 
is used. The potential energy per nucleon given by the 1st- 
, 2nd- and all-order ring diagrams are denoted respectively 
by PE-lst, PE-2nd and PE-ring. Eo - 1st and Eo - ring 
are given by the 1st- and all-order ring diagrams respectively. 
The experimental result is from [s^]. AH energies are in MeV. 



P/Po 


PE-lst 


PE-2nd 


PE-ring 


Bo-lst 


-Bo-ring 





-19.84 


-1.99 


-23.02 


-3.61 


-6.79 


0.50 


-20.67 


-2.28 


-23.60 


-4.44 


-7.37 


0.70 


-20.78 


-2.38 


-24.11 


-4.55 


-7.87 


0.75 


-20.79 


-2.40 


-24.25 


-4.56 


-8.01 


-Bo-exp 










-8.00 



E 

Q. 




may be noted that our results for p/po= 0.70 and 0.75 
are rather close to each other. We have calculated the 
potential enengy(PE) from the (^2^ + V^'"^) interac- 
tion. Results given by the 1st- and 2nd- and all-order 
ring diagrams, denoted by PE-lst, PE-2nd and PE-ring 
respectively, are listed in the Table. It may be noted that 
PE-2nd is much smaller than PE-lst. And the difference 
between (PE-lst -t-PE-2nd) and PE-ring is rather small, 
about 5%. Bogner et al. |34| have found that Viow-k 
is suitable for perturbative nuclear matter calculations. 
The good convergence shown in Table I suggests that this 
interaction is also suitable for perturbative calculations 
of finite nuclei. 

Having seen that the medium-dependent three-nucleon 
force V^"^ has given desirable results for nuclear matter 
and closed-shell nucleus ^^O, we shall now study its ef- 
fects on the nuclear effective interactions for the degen- 
erate sd one-shell and the non-degenerate sdpf two-shell 
model spaces. We use the same methods as described 
in [2^ for the calculation and application of these effec- 
tive interactions, except for the following difference: In 
[isf only the two-nucleon interaction V2N was employed. 
In the present work we include both V2N and Vg'JJf'^. We 
first calculate the low-momentum interactions from these 
interactions, and using them to obtain the Q-box dia- 
grams as shown in Figs. 1 and 2. The same V2N and 
^med interactions we used earlier in our calculations for 
nuclear matter and ^^O are employed. The degenerate 
sd and non-degenerate sdpf effective interactions are ob- 
tained from the Q-box using respectively the LS iteration 
method If 



20[ and the EKKO iteration method [21, 
The nondegenerate sdpf effective interactions can also be 
calculated using the KK iteration method [23h25,] . 

We first apply the above effective interactions with 
V^^"^ to nuclei ^^O and ^^F. Since V^^'^ is density de- 
pendent, we need to know the 'local density' felt by 
the valence nucleons in these nuclei. It should be small, 
but its precise value is rather uncertain. In the present 
work, we shall estimate p^ by comparing the density pro- 
file of ^^O and that for the valence nucleons. In Fig. 6 we 



+ 



FIG. 6: Density profile of O (solid line) and the radial 
distribution of the fourth power of the Od shell-model wave 
function (dashed line). See text for other explanations. 



4 

2 

L 

z 

z 

?-2 

z 

z 

" -4 
-6 



> 

CD 




-5 -4 



-3 -2-1 1 
Vnn [MeV] 



FIG. 7: Comparison of the sd-shell matrix elements of 
Veff calculated from V2N alone (NN) and from V2N + Vzn 
(NN+NNN). 



plot the density profile Pcoreix) of ^^O obtained with its 
wave function assumed to be a closed shell-model s*p^^ 
core. We also plot the distribution $(r) = ^0(f)^ where 
4>{r) is the radial harmonic oscillator wave function for 
^5/2, with the scaling parameter ^=0.1. As seen from 
the plots, the pair of c?5/2 nucleons reside primarily in 
the low-density region of Pcore- Depending on the av- 
eraging procedure employed, we have estimated p^ from 
Pcoreir) and 4>{r), obtaining values ranging from ^ 0.015 
to ~ 0.030/to~'^. We have considered another scheme to 
estimate py . The rms radius < r^ai > for the shell- model 
1^5/2 orbit is indicated by an arrow in the figure. A simple 
scheme to estimate py is to let it equal to Pcore{< fvai >), 
the value obtained in this way being ~ 0.025/m~'^. This 
value may be reduced if realistic s.p. wave functions are 
employed. The ^3/2 orbit is nearly unbound, and its rms 
radius should be considerably larger than that for d^/2. 
Assuming hoj —10 MeV for the ^3/2 orbit, its < ryai > 
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FIG. 9: Same as Fig. 8 except for ^**F. 



would be ^ 3.8/m giving ~ O.Ol/m^'^ for this orbit. 
We believe that a suitable range for py is from 0.1 to 
^ 0.2po (po = 0.16/to~'^). In the present work, we adopt 
Pv — O.15po7 to illustrae the effects of V^^^'^. 

In Fig. 7 the matrix elements of the sc?-shell effec- 
tive interactions calculated from V2N with and without 
V^'''^, denoted respectively as 'NN-f NNN' and 'NN', are 
compared. As seen, the two sets of matrix elements are 
rather similar to each other, with the magnitudes of those 
including V^''' being slightly larger. The low-lying spec- 
trum of ^^O and ^^F calculated with the above interac- 
tions are presented in Figs. 8 and 9. The main effect of 
the 'NNN' force is a small downward shift for the ground 
states of ^^O and ^^F, while leaving the other states 
largely unchanged. This trend is consistent with what 
we have observed for nuclear matter calculations (see Fig. 
4), namely the effect of Vg'JJf'* is slightly attractive at low 
densities. With the inclusion of thre 'NNN' force, the 
ground-state energy of ^^F is in good agreement with 



experiment while that for ^^O is slightly overbound. The 
agreements of the other states with experiments for ^^O 
are better than ^^F. 



Using the method outlined in section II and [25|, we 
have also calculated the effective interaction for a sdpf 
two-shell space using V2N with and without the inclusion 
of V^'^'^. In Figs. 8 and 9, the spectra of ^^O and ^^F ob- 
tained with the former effective interaction are also pre- 
sented ('two-shell', 'NN-HNNN'). As seen, the 'two-shell' 
and 'one-shell' results are rather similar to each other. 
To further study the effects of the three-nucleon force 
^iW^^ we have extended the calculations of Figs. 8 and 

9 to nuclei ^^O and ^^F, with results presented in Figs. 

10 and 11. Here, similar to ^^O and ^^F, the ground- 
state energies of ^^O and ^^F obtained with V2N only 
are both too high compared with experiments (see the 
'NN' columns). As shown by the 'NN-HNNN' columns, 
the inclusion of V^^'^ slightly lowers these energies. The 
'one-shell' and 'two-shell' results are rather close to each 
other for ^^O. But they are slightly different for ^^F, with 
the 'NN+NNN' ground-state energy noticebly lower than 
the 'NN' one. 

As shown in Figs. 8-11, although there are qualita- 
tive agreements between the calculated and experimental 
spectra of ^^O, ^^F, ^^O and ^^F, there are significant 
disagreements between them. For example, the orderings 
of the calculated low- lying states of ^^O are in fair agree- 
ment with the experimental ones, but their relative spac- 
ings are not. As shown by the figures, the effects from 
^sn"^ are generally small and attractive, mainly slightly 
lowering their ground-state energies while having nearly 
no infiuence on the relative spacings of the calculated 
spectrum. It should be useful to further study how to re- 
duce the above disagreements. We have employed a low- 
order approximation for the Q-box, including only 1st- 
and 2nd-order diagrams. The inclusion of higher-order 
diagrams r nay alter our results. As carried out by by 
Holt et al. [35|, certain classes of planar diagrams of the 
Q-box can be summed up to all orders using the Kirson- 
Babu-Brown induced interaction method. We plan to 
extend our present calculations by including these pla- 
nar diagrams, and study their effects to the shell-model 
description of these nuclei. ? 



IV. SUMMARY AND CONCLUSION 

We have carried out shell-model calculations for nuclei 
18q^ i^F, ^^O, ^^F and ^^O using effective interactions 
derived from the chiral two- and three-nucleon poten- 
tials. We start from the N^LO two-nucleon potential 
V2N together with an in-medium three-nucleon potential 
V^'^) which is a density-dependent effective interaction 
derived from the lowest-order chiral three-nucleon force 
V3N by Holt, Kaiser and Weise [13, Eil. The Viow-k 
low-momentum interactions are then derived from these 
potentials, and are used in a folded-diagram expansion 
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for calculating the V^jf interactions for both the degen- 
erate sd and non-degenerate sdpf model soaces. The 



well-known Lee-Suzuki iteration method |19l . |20[ is used 
for the derivation of the sd effective interactions, while 
the recently developed iteration method of Okamoto et 
al. 21| is employed for the sdpf ones. 

We have studied the effects of V^^'^, which is an ef- 
fective two-body interaction depending on the medium 
density p, at various densities near and below the nu- 
clear matter saturation density po- We first apply the 
VJ^'^'^ interaction to nuclear matter calculations for which 
the medium density is uniform and well defined. A ring- 
diagram method for nuclear matter has been employed. 
The low-energy constants Cd and Ce, of respectively 
the one-pion and contact terms of the leading-order chi- 
ral three- nucleon potential Vsn, are not well known. We 
have chosen these constants by requiring that their use 
in our nuclear matter calculations with V2N and V^'"^, 
which is derived from V3N, satisfactorily reproduces the 
saturation properties of nuclear matter. As indicated in 
Fig. 4, the inclusion of Vy/f^ has largely improved our 
nuclear matter results comparing with those given by V2N 



alone, and the resulting nuclear matter saturation prop- 
ertiesa are in good agreement with the empirical values. 

We then calculate the sd and sdpf effective interac- 
tions with the inclusion of the above V^"^, and apply 
them to nuclei ^^O, ^^F, ^^O and ^^F. An ambiguity here 
is the determination of the local density felt by the 
valence nucleons of these nuclei. This density should be 
small compared with the saturation density pQ of nuclear 
matter, but it is difficult to determine its value precisely. 
We have estimated this value based on the density profile 
of ^^O given by a simple s'^p^^ shell-model wave function. 
This is an approximation which should be further inves- 
tigated and improved upon. In our calculations for these 
nuclei, we have adopted the estimation of p^/po=0.15. 
As indicated by Figs. 8-11, the main effect of V^'^ is to 
lower the ground-state energies of these nuclei slightly, 
making them in better agreements with experiments than 
the energies given by V2N alone. The effects of V^'''' to 
these nuclei are general small, mainly because their va- 
lence nucleons are in a low-density medium. 

The above Vjj™^'^ has also been applied to a shell-model 
calculation of closed-shell nucleus ^^O. The nucleons in 
closed-shell nuclei are embedded in nuclear medium of 
much larger density than that for the valence nucleons 
mentioned above. For them the effect of V^^'^ is ex- 
pected be more important. Indeed this is confirmed by 
our calculations. As illustrated in Table I, the inclusion 
of this interaction has been essential for satisfactorily re- 
producing the empirical binding energy of ^^O. It should 
be useful and of interest to further study this in-medium 
three-nucleon potential by applying it to closed-shell nu- 
clei and those with a few valence holes. 
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